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We present a microscopic theory of zero-temperature order parameter and pseudospin stiffness re- 
duction due to quantum fluctuations in the ground state of double-layer quantum Hall ferromagnets. 
Collective excitations in this systems are properly described only when interactions in both direct 
and exchange particle-hole channels are included. We employ a functional integral approach which 
is able to account for both, and comment on its relation to diagrammatic perturbation theory. We 
also discuss its relation to Gaussian fluctuation approximations based on Hubbard-Stratonovich- 
transformation representations of interactions in ferromagnets and superconductors. We derive 
remarkably simple analytical expressions for the correlation energy, renormalized order parameter 
and renormalized pseudospin stiffness. 

I. INTRODUCTION 

Bilayer quantum Hall systems at Landau level filling factor v = 1 have broken symmetry ground states that can 
be regardecLeither as easy plane ferromagnets or as excitonic superfluids and have been extensively studied over the 
past decadaTLLj. Bilayer two-dimensional (2D) electron systems consist of a pair of 2D electron gases separated by 
a distance d (d ~ 100 A) which is comparable to the typical distance between electrons within each layer. When a 
strong magnetic field is applied perpendicular to the electron layers, the quantum Hall regime, in which macroscopic 
Landau level degeneracy quenches the kinetic energy and enhances the role of interactions in determining physical 
properties, is reached. Bilayer quantum Hall systems exhibit a rich variety of broken symmetry statesmi 2 ] depending 
on Landau level filling factor and on the relative importance of interlayer and intralayer Coulomb interactions. For 
filling factor v = 1, the ground state is completely spin polarized and spin excitations are gapped because of Zeeman 
coupling of the electron spins to the (strong) magnetic field, justifying the neglect of quantum spin dynamics. The 
only remaining degrees of freedom are the intra-Landau level cyclotron orbit centers and the discrete layer index which 
can be treated as a pseudospin label with "up" denoting a state localized in the upper layer and "down" denoting a 
state localized in the lower layer. At total filling factor v = 1 this system's ground state has spontaneous interlayer 
phase coherence (easy- plane ferromagnetism in the pseudospin language), and is incompressibletrLj. Since each layer 
has filling factor v — 1/2, a compressible-state filling factor for isolated layers, this quantum Hall effect is entirely 
due to correlations between the two layers. The system undergoes a phase transition at a critical layer separation d cr 
from this incompressible QHE state with_Dseudospin ferromagnetism to a disordered compressible stateoQ, possibly 
with other more exotic intervening statestil. 

In this paper we present a theory of quantum fluctuations about the Hartree-Fock mean-field-theory ordered ground 
state of the v — 1 bilayer. Our theory is based on an approximate expression for the ground state energy which includes 
quantum fluctuations and generalizes the random phase approximation to cases when both direct and exchange 
interactions are important. The ground state energy itself is of little direct physical relevance, however our approach 
is sufficiently flexible so that we can evaluate the dependence of the ground state energy on external fields and 
the winding wavevector of spiral states, to obtain quantities of physical interest. The most important macroscopic 
parameters characterizing order in quantum Hall bilayer pseudospin ferromagnets are the total pseudospin polarization 
and the pseudospin stiffness p, which characterizes the energetic cost of long wavelength fluctuations of the pseudospin 
field. In mean-field theory, all electrons occupy symmetric bilayer states, which corresponds to having all pseudospins 
polarized along the positive x-axis. We will use the pseudospin polarization normalized to its mean field value as a 
dimensionless order parameter, m x , which we will determine by evaluating the dependence of ground state energy on 
interlayer tunneling amplitude A t . The pseudospin stiffness will be determined by evaluating fluctuation corrections 
to the energy of spiral states, in which the pseudospin orientation varies spatially at a constant rate. 

At a formal level, the calculations described in this paper can be adopted to describe the influence of quantum 
fluctuations on the ground state of any easy-plane itinerant electron ferromagnet or any superconductor. In the case 
of a superconductor, the analog of the pseudospin stiffness is the superfluid density. The present work is motivated in 
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part by an interest in applications to these more general problems. The bilayer quantum Hall system is an ideal test 
case, both because its translational invariance permits many elements of the calculations to be performed analytically, 
and because the importance of quantum fluctuations in the ground state can be adjusted over a large range simply by 
adjusting the intcrlayer spacing. For zero layer separation, the bilayer systems is equivalent to a single-layer system 
with a spin degree of freedom and a Zeeman coupling A z that plays the same role as the tunneling amplitude At in 
spin-polarized bilayer systems. In the d — limit, the mean-field ground state is exactcl because the interaction term 
in the Hamiltonian is pseudospin invariant. At finite d, the difference between the intralayer Coulomb interaction, 
Va(q) = 2ire 2 /eq and the interlayer Coulomb interaction Ve(q) = VA(q)e~ qd reduces the interaction term's symmetry 
from SU(2) to U(l). Because of this reduced symmetry, only the z component of the total pseudospin polarization 
commutes with the interaction Hamiltonian. Hence, even though the mean-field-theory splitting between symmetric 
and antisymmetric single-particle state energies is enhanced by interactions, a fully pseudospin-polarized state cannot 
be an eigenstate of the microscopic Hamiltonian and in particular, cannot be the ground state. In the description we 
use, the ground state many-body wavefunction is a linear combination of states with collective excitations embedded 
in the fully-polarized Hartree-Fock mean-field ground state. 

The paper is organized as follows. In section [n| we describe the functional integral approach to this problem. It 
is knowu. that in these systems, both electrostatic and exchange fluctuations are essential to the collective mode 
physicsou, with the former dominating the cost of ^-direction pseudospin fluctuations and the latter controlling 
the easy-plane pseudospin stiffness. This situation contrasts with that of metallic ferromagnets in which exchange 
interactions, alone determine the collective mode behavior. No standard Hubbard Stratonovich (HS) transformation 
treatmentliia can deal with both simultaneously. Section [0] formally summarizes the generalized Hubbard-Stratonovich 
we employ that treats direct and exchange channels on an equal footing, and, establishes some of the notation we 
will use. Our approach is based on one developed previously by Kerman et alr3 Section HI formally summarizes the 
generalized random phase approximation (GRPA) for particle-hole correlation functions and its relationship to the 
generalized HS transformation. In section [V we apply our formalism to a double layer system at filling factor v = 1. 
We calculate the dispersion of collective modes, their contribution to the grand potential, and the effect of fluctuations 
on the pseudospin polarization and pseudospin stiffness. We conclude the paper with discussion in section |v|. 



II. AUXILIARY FIELD FUNCTIONAL INTEGRAL APPROACH 



In bilayer quantum Hall systems, due to the difference between interlayer and intralayer Coulomb interaction, there 
is a capacitive energy cost associated with charge-imbalance between the two layers. This anisotropy of the Coulomb 
interaction in pseudospin-space makes it necessary to treat the fluctuations in direct and exchange channels on .an 
equal footing. A general way to take into account such fluctuations is the Hubbard Stratonovich transformationliij. 
In this approach, we introduce an integral over a Bose field <f> which converts the fermionic two-body interaction 
in the original Hamiltonian into a one-body term coupled with the field </>. The Hamiltonian is then quadratic 
in the fermion operators so that the trace over these degrees of freedom can be evaluated exactly to obtain an 
effective action for the Bose field. Stationary phase approximations to the bosonic action yield mean-field theories. 
The nature of the mean-field state, for example whether it is characterized by a Hartree mean-field, an exchange 
mean-field, or a pairing mean-field depends upon the way the Hubbard Stratonovich transformation is done. The 
nature of the collective excitations that emerge when Gaussian fluctuations are allowed in the boson field also change 
qualitatively. For example fluctuations in the Hartree mean-field reveal plasmon collective modes, and fluctuations in 
the exchange mean-field of a ferromagnet reveal spin- wave collective modes. For a bilayer system, where the collective 
mode properties are determined by fluctuations in both channels, we need a method which treats both Hartree and 
exchange mean-fields on an equal footing. It is known that the standard HS transformation does not yield a Hartree 
Fock mean-field state and cannot capture both Hartree and the exchange fluctuations^. 

Kermanj— Levit, and Troudet have presented a generalization of the HS transformation which overcomes these 
limitationscj. To establish our notation (which differs from that of Kerman et at), we briefly review the previous 
work which describes an approach for obtaining systematic corrections to Hartree-Fock mean-field approximations 
for the grand potential. Our approach can also be_used to systematically improve upon Hartree-Fock mean-field 
approximations for the one-particle Green's functionlHi 

Consider a many-fermion Hamiltonian in second quantized form, 

a/3 ot(5~i5 

where if is a one-body (kinetic) term and V is a two-body interaction. Here a denotes (a set of) single-particle 
quantum numbers. For example, in the case of double layer systems a = (n, k, a) where n is the Landau level index, k 
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is the intra-Landau level index and a is the pseudospin index. Introducing pair-labels a = (a' a) and the density matrix 
p a = c^/Ca, the kinetic energy term becomes K a p a , where summation over repeated indices is assumed. Similarly, 
the interaction term can by written as V a b£l a b where fl a b =■ p a Pb '■ is the normal-ordered two body interaction and 
V a b — (a' P'\V\a(3) . Then in an obvious matrix notation, Eq. ( [2.l| ) becomes H = Kp + vCl/2. The partition function 
in the grand canonical ensemble is given by 



Tr 



aN-j3H 



lim Tr e 



nN 
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M 



eKp 3 - -VSlj 



(2.2) 



where T denotes the (imaginary) time-ordered product, N is the number operator, e = /3/M and we have retained 
only the imaginary time index. In the thermodynamic limit a/j3 — > p, the chemical potential. Since the number 
operator commutes with the Hamiltonian we can expand only the e~@ H term. The central idea of this generalized 
Hubbard transformation is to expand around an arbitrary two-body interaction U, thus treating Kp + Utt/2 as the 
dominant term and (V — U)Q/2 as a perturbation. Using the Gaussian identity 
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at each time index j, a typical term in Eq.( 2.2) becomes 
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Here Af 2 = AfxAf stands for a repeat sum over the N single-particle state labels in the Hilbert space and V(j)j implies 
both a product over single-particle labels and the relevant normalization factors. The term 4>jPj (Bose field coupling 
to the density operator) is odd in <j>j and does not contribute to the integral. The combination ^C/ -1 ^ / 'Af 2 gives 
a factor of 1/e after integration and we recover the interaction term. Thus, even though the interaction term has a 
prefactor e 2 and the kinetic term has a prefactor of e, both terms are of the same order. We stress that the limit 
e — > must be taken after all <j> integrals are done. Introducing these fields at each time-step j we get 
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(2.5) 
(2.6) 



Systematic approximations to the grand potential are obtained by expanding the bosonic action ( |2.6| ) around its 
minimum. In the limit as e — > 0, the configuration of fields 4>° which minimizes the action {dS / d4>\<f,a = 0) is given by 



4>%{j) = P««'^(fe'/3'(i))^, 

where we have defined the thermal average as 



,M; 



(2.7) 
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ehi 



Tre^rn^Cl-e/ii) 



(2.8) 



with the mean-field Hamiltonian hi = (K — <j)f)p- We recall that Eq.( ^J) is an exact identity and the interaction 
term with prefactor e 2 is of the same order as the kinetic term with the prefactor e after all the (^-integrals are done. 
Stationary phase approximation, however, corresponds to replacing the measure T>ipj by T>4>j5(4>j — ifi®). Hence in the 
stationary phase approximation, interaction term with the prefactor e 2 can be neglected and does not contribute to 
the mean-field solution. 

Henceforth let us assume a stati c mean-field solution so that the mean-field values 4>°j are independent of the time 
index j. It is clear from Eq.( 2.7) that the mean-field Hamiltonian is solely determined by the trial interaction U 
and is independent of the actual two-body interaction V. We get the Hartree (Fock) mean- field by choosing U = — V 
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(U = +V ex ), while — U = (V — V ex ) — V A (the antisymmetrized interaction) gives the Hartree-Fock mean-field. The 
mean-field grand potential is given by 



^0 = ^<S[<£°] = ~ o E fnfmUnmnm ~ ~g E ln(l + e° 



-far 







(2.9) 
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where h\n) = (K — (j) )p\n) — e n \n) describe the mean-field eigenfunctions and eigenenergies, and /„ = {l + e H 
are the Fermi occupation numbers. The single-particle representation is determined completely by the trial interaction 
U and is independent of the actual two-body interaction V. For the mean-field ground state energy we get 
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where the subscript h stands for occupied (hole) states. When U — — V, we get the direct contribution to the ground 
state energy, whereas only for U — —V A we recover the ground state energy as the expectation value of the original 
Hamiltonian. 

We improve upon the stationary phase approximation by considering quadratic fluctuations a roun d the mean field 
<f)° . These, of course, depend upon the true microscopic interaction V. Expanding the action ( 2J3) to second order 
in fluctuating Bose fields <j>j = <j)° + £j gives 
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where only time indices are explicitly shown. The matrices D and S are defined by 

D afaS = (P*l3ti)py6(3'))<l>° - (p a p(j)) <i> {Pys(j'))<p° 

and 
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(2-11) 



(2.12) 



(2.13) 



Here a, (3 are the single-particle labels and j, j' stand for the time indices. The fluctuation contribution to the grand 
potential is given by 
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(2.15) 



We have to be careful about taking the limit e — > since the dimension of the matrix M is of the order 1/e (see 
Kerman et al. for details). The resulting grand potential can be naturally divided into a quasiparticle contribution 
and a fluctuation contribution from the low-lying collective modes. The quasiparticle contribution is given by 
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We emphasize that the quasiparticle contribution reproduces the Hartree-Fock approximation for the ground state 
energy irrespective of the choice of trial potential U, even though the single-particle eigenstates are still determined 
by U. Among all the trial potentials, the antisymmetrized interaction U = —V A corresponding to the Hartree-Fock 
mean-field optimizes the quasiparticle grand potential, d^l QV ldV\ xj— -y A = 0. Hence Hartree-Fock mean-field is a 
particularly good choice for the trial potential. It follows from Eq. ( 2.11| ) that the fluctuation contribution is given by 
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(2.18) 
(2.19) 



where Ae mn = e m — e„ are the mean-field quasiparticle splittings, and w„ are the collective-mode energies obtained 
by diagonalizing the modified fluctuation matrix M — UM . Here, we have used that M is a symplectic matrix whose 
eigenvalues occur in pairs of opposite signs and 
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where iu k — 2Trk/(3 is a bosonic Matsubara frequency. The subscrip ts p, h in Eq.( p. 19 ) stand for unoccupied 
(particle) and occupied (hole) states, respectively. We can see from Eq.( 2.19] ) that the correlation energy is related 
to the differences between collective mode energies and the Hartree-Fock quasiparticle excitation energies. 

The functional integral formalism presented here is quite general and allows calculation of grand potential and 
collective modes around any mean field. Henceforth we will concentrate on expansion around the Hartree-Fock mean 
field which is accomplished by choosing U = ~V A . This choice gives, at the mean-field level, the exact ground state 
when the layer separation is zero. In the following section we discuss corresponding approximations in diagrammatic 
perturbation theory. 



III. GENERALIZED RANDOM PHASE APPROXIMATION 



In this section we describe the generalized random phase approximation (GRPA) for particle-hole response functions, 
which we relate later to the formal fluctuation energy expressions derived in the preceding section. Using a general 
notation, imaginary-time particle-hole response functions are defined by ^- abcd ( T ) = (Tc|(r)cb(r)cj(0)c t ;(0)). For non- 
interacting electrons these response functions have a simple formal expression in the representation of single-particle 
Hamiltonian eigenstates: 



n nm m i 
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(3.1) 



where we have Fourier transformed their imaginary time dependence and iu k = 2nk/f3 is a bosonic Matsubara 
frequency. These response functions have poles at the non- interacting electron particle-hole excitation energies, 
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In the generalized random phase approximation (GRPA), the single particle eigenstates are 



replaced with those obtained by diagonalizing the Hartree-Fock Hamiltonian hnF and collective fluctuations in direct 
and exchange potentials are captured by summing "ladders" and "bubbles" as summarized diagrammatically in Fig. ^. 
Including these corrections shifts the response function poles from differences of single-particle energies to collective 
excitation energies. Consider, for example, a typical ladder diagram having p interaction lines. In a convenient matrix 
notation, the value of such a diagram is given by 



(2p+i) terms 



X P (iu k ) = 'D(iu k )V ex D(iu k )---V ex D(iu k 
■b\V ex D\cd) = (bc\V ex \ad)D cd = (cb\V\ad)D cd 



(3.2) 
(3.3) 



Here the p th order diagram has been written in terms of the Hartree-Fock particle-hole response function D and the 
exchange interaction V ex . These diagrams capture the effect of exchange fluctuations around the Hartree-Fock mean- 
field. Adding these contributions for all p gives the ladder-sum approximation to the susceptibility. This geometric 
series can be formally summed to yield the following matrix for the ladder-sum response function xi 



X i(iu k )=D{iu k )-[l-V ex D{iu k )\ 



-l 



(3.4) 



Finally we consider a typical bubble diagram where instead of the Hartree-Fock response function D, the ladder-sum 
response function xi is used. In matrix notation, a term with m interaction lines becomes 
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X m {iu k ) = Xi( iu k)Vxi(iu k ) ■ ■ ■ Vxi(iu k ) ■ (3.5) 

" v ' 

(2m+i) terms 

Again, summing the resulting geometric series we get the following expression for the bubble-and-ladder sum suscep- 
tibility xu 

Xu^uk)- 1 = xiiiuk)- 1 +V = Diiuk)- 1 + V A . (3.6) 

In general, the determinant of the response function matrix xu has isolated poles and branch-cuts. By convention, we 
regard the well-defined collective modes due to the isolated zeros of det(l + V A D) as the collective modes, whereas 
zeros which lie in a continuum in the thermodynamic limit are associated with renormalized particle-hole excitation 
energies. Assuming that fluctuations here include only the collective modes, we get 



(3.7) 



det Xbi(iu k )= TT . 2 

-*■ - 1 - lui — CO 

where w v are the collective-mode frequencies determined by the (isolated) zeros of det(l + V A D). 

With these general expressions it is easy to see the equivalence between the functional integral approach and 
diagrammatics. Let us consider the collective-mode dispersions obtained from the response function (diagrammatics) 
and the quadratic fluctuations (functional integral approach). Using antisymmetrized interaction as the tunable 



potential, U = —V , the determinant of the modified fluctuation matrix M — UM [Eq.( 2.11)] becomes 



dct M = Y[ dct [l + V A D{iu k )] . (3.8) 

iuk 

Comparing ( |3~6[ ) and ( |3^) , we see that collective-mode energies obtained from these two methods are identical. In 
other words, quadratic fluctuations around a Hartree-Fock mean-field give the same result as the GRPA (bubbles and 
ladders) as far as collective-mode dispersion is concerned. If we had used an exchange-field (Hartree-field) standard 
Hubbard- Stratonovich transformation, the Gaussian fluctuations would have had resonances at the poles of xi (Xfc)- 
It is clear that whenever xu has a pole structure that is grossly different from that of both xi or Xb> neither standard 
HS approach will be adequate. Bilayer quantum Hall systems provide one example of such a circumstance. 

The relationship between the two approaches requires more care when approximating the grand potential. The 
venerable random phase approximation consists of summingpthe set of bubbles with n > 2 interaction lines and gives 
the following fluctuation contribution to the grand potentialHI 

n 6 = p Tr In [1 + VD] - \ tr (VD) . (3.9) 

This is the same as the grand potenti al ca lculated from quadratic fluctuations around the Hartree mean-field as can 
be checked by using U — — V in Eq.( 2.18). Similarly summing the "bubbles" with antisymmetrized interaction V A 
gives 

fl = — Tr In [l + V A D] - - tr (V A D) . (3.10) 

In the functional integral approach, this contribution is that obtained by considering quadratic fluctuations around 
the Hartree-Fock mean- field. The diagrammatic content of Eq.( 3. 10| ) is shown in Fig. |^. The set of diagrams on the 



left, (a), consists of vertex-corrected bubbles with n > 2 interaction lines. These diagrams can be formally summed 
and we obtain Trln [1 + Vxi] ~ tr(Vxj). This set of diagrams has been used-an the literature to approximate the 
grand-potential when both, direct and exchange, fluctuations are importantlii However the formal expression for 
the grand potential that follows from the functional integral approach includes another class of diagrams which are 
shown on the right in Fig. ^. The set of diagrams with n > 2 interaction lines, labelled (b), gives the contribution 
Trln [1 — V ex D] — tr(— V ex D), whereas summing the diagrams in the set labelled (c) we get tr(Vxi) — ti(VD). The 
equivalence of sum of these three contributions with the formal expression is most easily seen by usin g the fact 



that det [l + V A D] = det [1 + Vxi] • det [1 - V ex D]. In identifying the Feynman diagram content of Eq.( |3T0|) it is 
necessary to take care in using the pair-state matrix notation that simplified the formal developments in section O. 
As we will see, the inclusion of this set of diagrams substantially alters the renormalization of the system parameters. 
We refer to the sum all three classes of diagrams as the GRPA. 
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IV. APPLICATION TO DOUBLE LAYER SYSTEMS 



In this section, we apply the formalism developed above to bilayer quantum Hall systems at total filling factor 
v = 1. We neglect spin so that the pseudospins and the orbit centers within a Landau level are the only dynamical 
degrees of freedom. Due to the properties of Landau level wavefunctions, there existe-a unitary transformation which 
allows us to associate a two-dimensional wavevector with a pair of orbit center labekslS. In other words, momentum is 
a good quantum number for the particle-hole pairs. This property is responsible for the considerable progress that can 
be achieved analytically in the following calculations. We want to calculate the correction to the mean-field ground 
state energy (or grand potential) due to collective excitations for spiral states and ferromagnetic states with interlayer 
coupling. The renormalized order parameter m x and the pseudospin stiffness p are obtained from the ground state 
energy by taking appropriate derivatives. 

The Hamiltonian for a double layer system with interlayer tunneling amplitude At is 

k,<7 

V = - ^{^^2^^ + a 1 a 2 V x \k 3 a 1 k i a 2 )c\ iai c\ 2a2 c k ^ 2 c k3< y 1 . (4.2) 

Here, a =|= — ~o denotes a state in top layer, while a =\ denotes a state in the bottom layer. The ki the are 
y-components of the canonical momenta which are good quantum numbers of the single-particle Hamiltonian in the 
Landau gauge A = (0, Bx,0), V = (Va + Ve)/2 and V x = (Va — Vg)/2 are sums and differences of interlayer and 
intralayer Coulomb interactions. V x characterizes the anisotropy of the interaction in the pseudospin space. The 
energetic splitting between symmetric and the antisymmetric Hartree-Fock eigenstates is enhanced by interactions, 
A s as = A t + A sb , where 

A sb = T E (0), (4.3) 

T ^ = J E ^(Ple-^V-^ 2 (4.4) 
p 

and A = 0,x, A, E. The T\(q) are interactions between an electron and an exchange holeEl whose center is separated 
by a distance ql 2 . 

First we concentrate on the case with zero interlayer tunneling. If the external tunneling is absent, A t = 0, there is 
a rotational symmetry in the x — y plane in the pscudospin-space. In this case the Hartree-Fock mean-field equations 
support a class of "tumbling" or spiral order-parameter solutions m(x) — (cosQx,sinQx). Ina superfluid language, 
this corresponds to a supercurrent solution parameterized by a pairing wave-vector Q = Qixli. There are a series of 
metastable spiral ground states with different values of Q. The Q-dependent eigenstates and eigenenergies are given 
by 

|fc,±} = 4jlT)±e 4QfeZ2 |i)l, (4.5) 







V2 

r B (Q)/2 = Aq/2 = -e° + . (4.6) 



For Q = we get the symmetric and the antisymmetric states as eigenstates and the spiral state reduces to the 
ferromagnetic state. By evaluating the fluctuation correction to the ground state energy at all values of Q we will 
be able to estimate fluctuation corrections to the spin-stiffness of bilayer quantum Hall fcrromagnets. The mean-field 
Green's function is diagonal in the Hartree-Fock eigenstate basis. The Hartree-Fock susceptibility Xhf ^ = is 
therefore given by 

where Sn = {n + — n_) is the difference between Fermi factors. As expected the Hartree-Fock susceptibility diverges 
at the quasiparticle energy gap Aq. Due to the absence of dispersion in Landau band energies, the ladder-sum for 
the susceptibility can be evaluated analytically. We find that 
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Here we have introduced 
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(4.9) 



r ++ (q) = [2T A (q) + T E (q- 
T + _(q) = [2T A (q)-F E (q- 
£ 2 = (A Q - r ++ ) 2 - 



• Qx) + T E (q - Qx)} /4, 

• Qx) - Y E {q - Qx)] /4, 



(4.10) 
(4.11) 
(4.12) 



Note that the ladder-sum susceptibility \i nas a P°l e a t i^n = £{q)- It is clear from Eq.( |4. 1 2| ) that this pole arises due 
to exchange fluctuations. In the functional integral approach, iVL n = £(q) denotes the zeros of det (1 — V ex D). Since 
all other elements of the susceptibility matrix vanish, wc will henceforth use only two indices to denote the matrix 



elements, x 



+ h = 



X 



= \ ++ an( I so on. Furthermore, since x + (?)*^) = X + (<f; — *^) an d X ++ = X 



X 

there are only two independent matrix elements. 

We include the bubble diagrams to take into account electrostatic effects. This gives the following equation for the 
susceptibility matrix Xbi(q, i^n) 



Xu 
Xu 



Xi 

xt 



+- 



+ (-2nl 2 e^ 2 / 2 v^) ■ 



(xr 

(xt + 



xT 

Xi 



(xt 
(xt 



xT 

Xi 



Xu 
Xu 



(4.13) 



Here 2-Kl 2 v^(q) = e i 2 * 2 1 2 [Va(<t) — V E (q) cos Qq y l 2 ] /2 represents the electrostatic fluctuation effects. This GRPA 
susceptibility has a pole at the the collective-mode energy given byO 



E sw(q,Q) = a Q(q)- b Q(q), 

a Q (q) = A Q -r A (q)+2v$(q), 

b Q (q) = A Q - X - [T E (q + Qx) + T E (q- Qx) 



(4.14) 
(4.15) 

(4.16) 



In the functional integral approach the same result for collective mode dispersion is obtained from the zeros of 
determinant of the fluctuation matrix M = (l + V A D) . Since momentum and frequency are good quantum numbers, 
the fluctuation matrix is diagonal in these indices, and effectively only has pseudospin labels. In the eigenstate 
representation, the interaction matrix elements are given by 



(k 1 a 1 ,k 2 (T2\V\k 3 <J 3 ,k4<T 4 ) = i(fcifc 2 |K4|^3^4) [1 + 0-102 (730-4] 



(4.17) 



+ -7{kik 2 \V E \k 3 k4) 



aicr 3 e 



iQ(k 3 -k!)l 2 



a 2 a i e 



iQ(ki-k 2 )l 2 



where fcj are the intra-Landau level indices and o-j = ±1 denote the fc-dependent pseudospin eigenstates. Further 
constraints on the pseudospin labels are imposed by the fact that the Hartree-Fock susceptibility is nonzero only when 
the particle-hole pseudospins are opposite. Hence, the relevant part of the fluctuation matrix becomes 



(a 1 a 2 \l + V A (q)D(m n )\a 3 ,a i ) 



l + (v§- r ++ ) 



$11 



Sn 



-); 



6 11 



1 



Sn 



(4.18) 



We emphasize that this 2x2 submatrix of the full 4x4 fluctuation matrix has the same determinant as the entire 4x4 
matrix. Hence it is sufficient to evaluate this submatrix as far as the contribution to the grand potential is concerned. 
In other problems, for example, when considering the effect of collective modes on the one-particle Green's functionES, 
the entire collective-mode propagator is required. Notice that in the functional integral approach, the effect of Hartree 
fluctuations appears via the use of antisymmetrized interaction U = —V A instead of exchange interaction U = V ex . 
Since a spiral state spontaneously breaks the rotational symmetry in the x — y plane in pseudospin space, collective 
modes are pseudospin- waves where the pseudospin polarization deviates slowly from the spiral-state polarization. The 
term aq(<z) in Eq.( 4.15) expresses the energy cost of fluctuations out of the x — y plane. These fluctuations correspond 
to moving charge from one layer to the other and are finite even at zero wavevector, a,Q(q — 0) ^ 0. In contrast the 
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term 6q(<?) represents the cost of pseudospin rotations in the x — y plane. In particular the fact that bq(q = 0) = 
indicates the Goldstone nature of these collective excitations. We stress that the preceding results for spiral-state 
susceptibility matrix and collective-mode dispersion are valid only when there is no external tunneling between the 
layers. 

Figure ^ shows a typical plot of collective- mode energy E sw (q,Q). We emphasize that the anisotropic nature of 
the dispersion is due to the choice of the spiral wavevector Q = Qx and the non-monotonic nature of the energy is 
due to the competition between the Hartree and exchange terms in aq{q). As q — > o o the collective- mode energy 



4.14) reduces to the results 



approaches the quasiparticle gap Aq_ For the uniform case (Q = 0) -the dispersion 
previously obtained by diagrammaticaj and single-mode approximation!! 

Now let us consider the explicit expressions for contributions to the grand potential due to various sets of diagrams 
mentioned in the previous section. The functional integral approach which we have discussed in previous sections 
gives the following fluctuation contribution to the grand potential 



a. 



— In TT 

2/3 11 



in 2 n -E 2 sw (p,Q) 



ml 



A 2 

Q 



- \ tr(V A D) 



(4.19) 



where the first term - det (l + V A D) - has been expressed in terms of the collective- mode frequencies E sw and the 
quasiparticle gap Aq [see Eq.( 2.2C| )]. In contrast, the contribution to the grand potential obtained by summing 
bubbles with vertex corrections - the approximation which has been used in the literatureEj - is given by 



in 2 n -E 2 sw (p,Q) 
iW n -£ 2 (p,Q) 



- 2 M^xO 



(4.20) 



It is clear from Eqs.( 4.19 ) and ( 4.2C ) that the two approximations are different. 

We now calculate the effect of these collective modes on pseudospin stiffness and pseudospin polarization of the 
system. These parameters are obtained from the ground state energy Eq{Q) by taking appropriate derivatives 



P = 



1 d 2 E G (Q) 
A dQ 2 1 
AttI 2 dE G {0) 



=0, 



(4.21) 
(4.22) 



In the following we use Eqs.( 2.17 ) and ( 2.19| ) to approximate the ground state energy. For the spiral-state case, 
neglecting the kinetic contribution (which is Q-independent), we obtain 



E, 



1^(0) AQ 

2 2 



A 4tt/ 2 



4ttZ 2 



(4.23) 
(4.24) 



The Hartree-Fock stiffness is obtained from Eq.( 4.21) when the ground state energy is approximated by E qp . The 
renormalized stiffness is obtained b y usi ng Eq = E qp + E c . Note that the only Q-dependence in Eq is through the 
term [E sw (p, Q) — Aq]. Using Eq.( 4.21) we obtain the following explicit expression for the renormalized pseudospin 
stiffness, including correlations 



-y 

•A ^ 



2A 



a'pb p 



a p b' p ' 



2E sw (p) 



8ttI p HF 



(4.25) 



This equation is one of the principle results of this work. Here a p = aQ=o(j>) and so on, the primes denote second 
derivative with respect to the wavevector Q and we have used the fact that Aq =0 = —8tt1 2 phf- Figure ^ compares 
the Hartree-Fock and the renormalized pseudospin stiffness. At d — the Hartree-Fock ground state is exact and as 
expected, there is no stiffness renormalization. We note that for small layer separations d < I fluctuations enhance 
the pseudospin stiffness, a possibility that has not been anticipated previously. In our theory a p vanishes at a finite 
value of p > Z _1 at layer separations that exceed a critical value d cr . In Hartree-Fock theory—the ground state changes 
from a uniform coherent state, to a coherent pseudospin-density wave state at this pointlij. This instability point 
has previously been identifiecfl, on the basis of heuristic arguments, with a phase transition at which pseudospin 
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ferromagnetism is lostE2. As a p approaches zero, the first term in Eq.( 4.25) becomes large and negative and p reaches 
zero for d slightly smaller than d cr . While our Gaussian fluctuation calculations are certainly not systematic in this 
regime, our result that the stiffness goes to zero rapidly for d near d cr is consistent with numerical exact-diagcpalization 
estimatesEl and recent suggestions that pseudospin ferromagnetism could be lost via a first order transitioned. 

To calculate the order parameter renormalization we start with a system having a finite tunneling amplitude A t . 
In this case the pseudospin symmetry in the x — y plane is explicitly broken, and the Hartree-Fock mean-field state 
is uniformly polarized along the ir-axis. The quasiparticle gap Asas includes the external tunneling contribution At 
and the collective mode is gapped at zero wavevector. Using the finite tunneling and zero spiral wavevector version 
of Eqs.( 



4.23), ( 4.24) we get the following expression for the renormalized order parameter 



2nl 2 



E 



(4.26) 



where e p = (a p + b p )/2. This anal ytic expression is the second principle result of this calculation. The two terms 



on the right hand side of Eq. ( |4.26| ) have a simple interpretation. First term is the Hartree-Fock result for the order 
parameter (m x — 1) and the second term represents its suppression due to collective excitations. Figure |^ compares 
the Hartree-Fock and renormalized order parameter. When d = 0, we have e p = E p and m x — 1. As d — > d cr , collective 
modes mix more strongly into the mean-field ground state and the order parameter is suppressed. We mmark that 
these renormalized parameters are qualitatively different from the parameters obtained in the literaturet£l by using 
Eq.( 4.20) as an approximation for fluctuation correction to the grand potential. 



V. DISCUSSION 



We have presented a functional integral approach to systematically calculate fluctuation corrections to a mean-field 
approximation for the grand potential. These corrections represent the contribution due to collective modes around 
an ordered mean-field ground state. Our approach takes into account fluctuations in both, direct and e xchange, 
particle- hole channels. The same expression for fluctuation corrections to the grand potential, Eq.( 2.18 ), can be 
systematically derived by coupling constant integration from GRPA particle-hole response functions. Remarkably, for 
bilayer quantum Hall systems it is possible to calculate these corrections analytically. The pseudospin stiffness and 
order parameter for the phase-coherent state are obtained from grou nd st ate energ y by t aking appropriate derivatives. 
The two central results of our work are the explicit expressions, Eq.( 4.2£ ) and Eq.( 4.26), for renormalized pseudospin 
stiffness p and the renormalized order parameter m x respectively. We note that these expressions, derived from a 
microscopic treatment, are very similar to those obtained by doing a linearized spin-wave analysis of an easy-plane 
XY spin model on a lattice. In case of a lattice model, though, the two-dimensional momentum sum would be cut 
off at p = 2t: i 'a where a is the lattice constant. 

Our approach presented here is quite generic and applicable to systems where fluctuations in more than one channel 
dominate the low-energy physics. A relevant example is the physics of underdoped cuprates where, presumably, 
fluctuations in the superconducting and antiferromagnetic channels are important to low-energy physics. The approach 
presented here can be used to calculate, for example, the renormalization of superfluid density or penetration depth 
due to quantum fluctuations. 

In this paper, we have presented spiral states as solutions to the mean-field equations when the interlayer tunneling 
is absent, A t = 0. In such states, the relative phase between single-particle states localized in the top and the bottom 
layer varies as e % Q x . These spiral states are also realized as mean-field ground states in bilayer systems with nonzero 
interlayer tunneling in the presence of an in-plane field B\\p. In the presence of such a field, the Aharonov-Bohm 
phase picked up by an electron tunneling from top layer to the bottom layer varies precisely in the same manner with 
wavevector Q = dB\\/l 2 B. 

The results presented here can be easily generalized to include the form-factors arising from finite well-widths. It 
is also possible to take into account the effects of finite temperature by using grand potential instead of the ground 
state energy to calculate the renormalized parameters, though it will require some more work. 

In this work we have limited ourselves to Gaussian fluctuations around the Hartree-Fock state. This means we 
are treating the collective modes as non-interacting bosons. If we were to go beyond the Gaussian approximation 
and include higher-order fluctuations, they would generate interactions between these collective modes, which are 
certainly important as the layer separation approaches critical layer separation, d — > d cr . A self-consistent treatment 
of fluctuations is required in this regime. Since we know self-consistent approximation schemes for spin-models on 
a lattice, it is an interesting problem to develop effective spin-models which incorporate the microscopic physics of 
bilayer systems. In the present treatment, the strong renormalization of the macroscopic parameters only close to 
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the phase-boundary indicates that the Gaussian fluctuation approximation is a reasonable approximation over a large 
range of interlayer separation. 
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FIG. 1. GRPA for the susceptibility matrix. The ladder-sum susceptibility \i takes into account exchange-potential fluctu- 
ations. The electrostatic fluctuations which are important in bilayer systems are taken into account by summing bubbles, Xbl 
is the GRPA susceptibility. We have written the ladder sum as a bubble sum with an altered interaction matrix element, in 
keeping with the algebraic structure of our formal development in a matrix representation with pair excitation labels. 
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FIG. 2. Diagrammatic content of Eq.( 3.1C ). The set of diagrams on the left, (a), is a random-phase-approximation-like 
bubble sum, but here the bubbles have ladder diagram vertex corrections. Only diagrams with n > 2 bubbles appear in this 
set. The symmetry factors of these diagrams match those obtained by expanding the In as usual. The sets of diagrams on the 
right, labelled (b) and (c), appear naturally in the functional integral formalism. All these sets of diagrams can be summed 
explicitly for double layer quantum Hall systems. 
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FIG. 3. Typical collective mode dispersion for the spiral state. The anisotropy in the dispersion is related to the choice of 
the spiral wavevector Q = Qx. The non-monotonic behavior is due to competing Hartree and exchange fluctuations present 
in this system. In our functional integral approach, the spin- stiffness can be related to the dependence of the collective mode 
energies on the spiral wavevector. 
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FIG. 4. Renormalization of the pseudospin stiffness due to quantum fluctuations. The stiffness is enhanced by fluctuations 
at typical values of d. The stiffness vanishes rapidly close to the phase boundary d cr as discussed in the text. Note the absence 
of renormalization at d = 0, when the Hartree-Fock state is the exact ground state. 
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FIG. 5. Dependence of the dimensionless order parameter on layer separation. The mean-field polarization is independent 
of layer separation d and is not sensitive to the changes in collective mode energies that occur as d approaches d cr . The 
renormalized polarization vanishes rapidly close to the phase boundary d cr . 
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